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We perform a detailed analysis of the importance of causality within the strong-field approxi- 
mation and the steepest descent framework for the recoUision-excitation with subsequent tunneling 
ionization (RESI) pathway in laser-induced nonsequential double ionization (NSDI). In this time- 
delayed pathway, an electron returns to its parent ion, and, by recoUiding with the core, gives part 
of its kinetic energy to excite a second electron at a time t' . The second electron then reaches the 
continuum at a later time t by tunneling ionization. We show that, if t' and t are complex, the 
condition that recollision of the first electron occurs before tunnel ionization of the second electron 
translates into boundary conditions for the steepest-descent contours, and thus puts constraints on 
the saddles to be taken when computing the RESI transition amplitudes. We also show that this 
generalized causality condition has a dramatic effect in the shapes of the RESI electron momentum 
distributions for few-cycle laser pulses. Physically, causality determines how the dominant sets of 
orbits an electron returning to its parent ion can be combined with the dominant orbits of a second 
electron tunneling from an excited state. All features encountered are analyzed in terms of such 
orbits, and their quantum interference. 
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I. INTRODUCTION 

Phenomena occurring in the interaction of matter with 
intense, low-frequency laser fields, such as high-order 
harmonic generation (HHG), above-threshold ionization 
(ATI) , and laser-induced nonsequential double ionization 
(NSDI), owe their existence to laser- induced recombina- 
tion or rescattering processes, in which a raeviously re- 
leased electron interacts with its parent ion [l| . The iden- 
tification of their common dynamical origin has been in- 
strumental to exploit these effects in a number of applica- 
tions, such as for the generation of attosecond pulses [2|| 
and the dynamic imaging of matter with subfemtosecond 
precision [3] . Since this mechanism is intimately linked to 
specific dynamical pathways, it also led to efficient semi- 
analytical approaches, such as the strong-field approxi- 
mation (SFA) (for seminal work on ionization and laser 
induced rescattering see, e.g., Q and [5[, respectively). 
In this approach, the quantum-mechanical transition am- 
plitude is associated to the orbits of an electron returning 
to its parent ion. The transition amplitude correspond- 
ing to a particular strong-field process takes the form of 
a multiple integral, which, in many cases, can be solved 
employing saddle-point methods !6|. 

It is also a well-known fact that the orbits along which 
the active electron returns typically occur in pairs. In 
particular, thresholds can often be associated to condi- 
tions where the two orbits of a given pair become almost 
degenerate. This violates the saddle-point assumption 
which leads to expressions in terms of individual orbits, 
but can be treated successfully in a uniform approxi- 
mation which describes the two orbits collectively tVi]. 
This approximation has been first applied in strong-field 
physics to ATI Q and, since then, has been used in a 
wide range of phenomena, such as HHG 9] and NSDI 



|10l - ll2| . This latter phenomenon is a typical example 
of a highly correlated two-electron process occurring in 
strong laser fields. The physical mechanism behind it is 
a three-step process in which the first electron is released 
in the continuum by tunneling or multiphoton ionization, 
and gains kinetic energy from the field. Subsequently, it 
is driven back by the field towards its parent ion, with 
which it rescatters. In this recollision, part of its kinetic 
energy is transferred to the core, so that a second electron 
is freed. 

There are several pathways through which NSDI may 
occur. Both electrons may, for instance, be released si- 
multaneously in a scattering process in which the first 
electron, upon return, provides the second electron with 
enough energy for it to overcome the second ionization 
potential. This process is known as electron-impact ion- 
ization. This pathway has been successfully used to ex- 
plain a number of experimentally observed features in 
the electron momentum distribution, such as peaks at 
non- vanishing momenta |lll . Il3l . Il4| , and the recently ob- 
served V-shaped structure which can been associated to 
the long-range electron-electron repulsion [13, [l3 ■ 

Apart from that, the second electron may be released 
in a time delayed pathway, in which the first electron 
promotes the second electron to an excited bound state. 
Near the subsequent field maximum, the second elec- 
tron tunnels from this excited state. These pathways are 
becoming increasingly important for two main reasons, 
which are directly related to attosecond-imaging appli- 
cations. First, the complexity of studied NSDI targets 
is systematically increasing, and with this internal exci- 
tations become more important. Second, delayed path- 
ways govern the below-threshold regime, for which the 
driving-field intensities are too low for the second elec- 
tron to leave by direct ionization |16| . 



These time-delayed pathways are far less understood 
[17| . Besides an overall controversy rela ting to the physi- 
cal mechanisms behind them (see, e.g., [la] for a detailed 
discussion of this topic), they can pose a fundamental 
problem associated to causality, which we address in the 
present work. A particularly relevant time-delayed path- 
way in NSDI is the "recollision with subsequent tunnel- 
ing ionization" (RESI). In previous work 19] it has been 
shown that this pathway can be understood as a rescat- 
tered ATI-like process for the first electron, followed by 
direct ATI for the second electron [23] • RESI is a non- 
standard case within strong-field physics, as the uniform 
approximation mentioned above requires some modifica- 
tions. In fact, a rigorous treatment of RESI within a 
saddle-point framework is a non-trivial problem as the 
rescattering of the first electron must precede the ioniza- 
tion time of the second electron. So far, this issue has 
been analyzed mainly in a classical framework, for which 
both times are real and thus can be ordered. In the SEA, 
however, the associated orbits have complex ionization 
times, which are required to account for the non-classical 
effect of tunneling. 

In this work, we address this complication of casual- 
ity for complex times and investigate how it affects the 
momentum distributions in the RESI pathway of NSDI. 
We approach the problem from the perspective of asymp- 
totic expansions, and resolve it by explicit construction of 
steepest-descent integration contours. The consequences 
of causality become apparent when one contrasts the case 
of a purely monochromatic (and hence infinitely long) 
driving field to the more realistic scenario of few-cycle 
driving pulses. 

This article is organized as follows. In Sec.|ll]we briefly 
recall the expression for the RESI transition amplitude 
and discuss the saddle point trajectories with complex 
rescattering and ionization times. In Sec. IIIII we calcu- 
late the individual momentum distributions of the first 
and the second electron for monochromatic driving or a 
few-cycle pulse, disregarding the correlations imposed by 
causality. In Sec. IIVI we describe how the causality re- 
quirement refiects itself in the complex time plane. This 
description is then used in Sec.|V]to compute correlated 
two-electron momentum distributions. Section I VI I con- 
tains our conclusions. 



II. BACKGROUND 

A. RESI transition amplitude 

The strong-field approximation (SEA) transition am- 
plitude for RESI with final electron momenta pi and P2 
reads 0,[2l|,[^. 
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Physically, Eq. ([T|) is associated to a rescattering pro- 
cess in which an electron, initially in a bound state JV'i ) 
with energy Ej^ , tunnels at time t" into a Volkov state 
|k(t")). Erom the time t" to the time t' , this electron 
propagates in the continuum, until it is driven back to its 
parent ion. Upon return, the electron scatters inelasti- 
cally with the core at time t' and, through the interaction 
Vi2, elevates the second electron from the ground state 
IV'2 ) of the singly ionized species (with energy E2 ) to 
the excited state \tp2 ) (with energy E^ )■ Einally, at 
a later time t, the second electron is released by tunnel- 
ing ionization from the excited state jV'e ) into a Volkov 
state |p2 (t)). Here, k(t") = k and P2{t) — P2 h^ the 
velocity gauge, and k(i") = k + A(t"), p2{t) = p + A{t) 
in the length gauge [281]. All the information about the 
binding potential V(ri) of the first electron and Mon(i"2) 
of the second electron, and the interaction Vi2(ri,r2) of 
the first electron with the core, are embedded in the form 
factors ([3]) , fl]) and ^ respectively, which we will assume 
to be constant over the parameter range in question. 

Throughout this work, we will consider both a 
monochromatic, linearly polarized field, for which 



A{t) = 2y^sm{ujt)i^ 



(6) 



with ponderomotive potential Up, and a few-cycle pulse, 
for which 

A(i) = 2v/C^sin2 f^j sin((^ + oot)e,, (7) 



where e^ denotes the polarization vector, N the num- 
ber of cycles in the pulse and the carrier-envelope 
phase. The associated electric field is defined by E(t) = 
— 9t A(t). A monochromatic wave is a reasonable approx- 
imation for long pulses [12| , or for few-cycle pulses when 
the carrier-envelope phase is integrated over. 



B. Saddle-point equations 

For large driving-field intensities, Eq. ([T]) is a 
strongly oscillatory integral which can be evaluated using 
steepest-descent methods [20] • This requires to obtain 
the saddle points where the action ([2]) is stationary. 



to the core in order to excite the second electron. One 
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and establishes a direct link to semiclassical orbits with 
complex times and actions. These orbits can be obtained 
efficiently by recognizing that the action splits into two 
independent parts, 

5i(pi, k, t", t') = Ei'h" + (4"^ - E!f^)t' 



for the first electron and 
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for the second electron. 



1. First electron 

Explicitly, the stationary conditions upon 5*1 lead to 
the equations 
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Condition (fTTj) states the conservation law of energy 
for the first electron when it reaches the continuum by 
tunneling. Condition (J12p constrains the intermediate 
momentum of the first electron so that it returns to the 
site of its release, and also guarantees that the inter- 
mediate momentum k of the electron is parallel to the 
laser field. Condition ([T3| gives the conservation of en- 
ergy upon rescattering of the first electron, and states 
that the final kinetic energy of the first electron is its ki- 
netic energy upon return, minus the energy it transferred 
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should note that, if E^ = E^ , the elastic rescattering 
condition for high-order above-threshold ionization is re- 
covered. 

In particular, the solutions of the saddle point equa- 
tions generally lead to complex times, as Eq. ([Tl]) admits 
no real solutions. This is a consequence of the fact that 
tunneling is not a classically allowed process. The imag- 
inary part of t" will be directly related to the width of 
the barrier through which the electron tunnels: the nar- 
rower the barrier, the smaller Im[i"] and the larger the 
tunneling probability. Furthermore, for given values of 
momentum, rescattering may or may not have a classi- 
cal counterpart, depending on whether the maximal ki- 
netic energy of the first electron upon return is larger or 



(8) smaller than the excitation energy ii^exc = £-2 ~ -^' 



p(e) 



In 

order to analyze this aspect, it is useful to recast Eq. p^ 
in terms of the electron momentum components piu ,pii_ 
parallel and perpendicular to the laser-field polarization. 
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This implies that a nonvanishing perpendicular momen- 
tum pij_ shifts the energy the first electron must provide 
to the core in order to excite the second electron. The 
maximal kinetic energy upon return has to be larger than 

an effective excitation energy E^cxc = (£-2 ^^2 )+-Pi_l/2 
for the rescattering process to have a classical counter- 
part. For that reason, one expects that the classically al- 
lowed region in momentum space will be most extensive 
for vanishing perpendicular momentum pi± [19l. |21|. As 
a function of the parallel component of momentum, the 
most favorable rescattering conditions are then achieved 
when the first electron leaves with piu — —A{t') around 
a maximum of |A(i')p, of which there are two per field 
cycle, and returns near the subsequent field crossing. 

In terms of momentum constraints, this translates into 
the condition {pi\\,Pi±) = {—A{t'),0) for the first elec- 
tron around which the partial momentum-space maps in 
the {pi\\,Pi±) plane discussed in Sec. IIIII will be centered. 
For the first electron, the solutions of the saddle-point 
equations will occur in pairs. These pairs correspond to 
the "long" and the "short" orbit of an electron rescatter- 
ing with its parent ion [23] . Each cycle will then contain 
two pairs of orbits, i.e., four orbits altogether. If a clas- 
sically allowed region is present, the rescattering times t' 
for each pair have vanishingly small imaginary parts of 
opposite signs, i.e., they are located either in the lower or 
upper complex half plane related to t' . The start times t" 
will always exhibit nonvanishing and positive imaginary 
parts. As one moves away from the center of a classically 
allowed region, the saddles in a pair approach each other 
closely, until they almost coalesce at its boundary. If, 
on the other hand, the parameter range is such that no 
classically allowed region is present, Im[t'] will no longer 
be vanishingly small in any momentum region. In this 
case, the physically relevant saddle will be located in the 
upper half plane. The remaining saddle will lead to ex- 
ponentially increasing results and must be discarded (for 



a detailed discussion see [21|). Longer pairs of orbits, in 
which the electron returns after having spent over a cycle 
in the continuum, will contribute much less to the yield 
due to the spreading of the electronic wave packet and 
therefore will be ignored. 



2. Second electron 
The stationarity condition upon S2 (p2 , t) yields 



[V2+M^? = -'^E^2\ 



(15) 



which, physically, gives the energy conservation upon ion- 
ization of the second electron. This final ionization pro- 
cess is classically forbidden throughout, as it always in- 
volves tunneling. If written in terms of the parallel and 
perpendicular electron momentum components p2\\, P2±, 
Eq. (Uni) reads 



[P2\\+A{t)f = -2E^^-Pl 



(16) 



This implies that a non-vanishing transverse momen- 
tum P2X effectively widens the potential barrier through 
which the second electron must tunnel. A direct con- 
sequence will be an overall decrease in the yield with 
increasing transverse momentum. The electron tunnels 
most probably at the field peak, and with p2± = 0. 
Hence, the momentum-space conditions for which ion- 
ization of the second electron is most probable read 
{p2\\iP2i) = (0,0). For the second electron, there ex- 
ist two saddles per field cycle, which do not coalesce. For 
vanishing momenta, the real parts of the corresponding 
ionization times lie at subsequent maxima of the laser 
field, half a cycle apart from each other. As the parallel 
electron momentum increases these saddles move towards 
the field crossings. Since, however, their contributions 
occupy the same region in momentum space, quantum 
mechanically they interfere. 

Equation ([15]) is identical to that governing direct ATI. 
The second electron will therefore obey similar momen- 
tum constraints as in this process (for which, however, 
causality does not play a role). For instance, for a 
monochromatic field, momentum-resolved distributions 
should be bounded by |p2|| | ^ 2-y/C7p, which corresponds 
to the traditional ATI cutoff energy of 2Up |19| . 



C. Causality and partial transition amplitudes 

The saddle point analysis described in the previous sec- 
tion provides two independent sets of orbits, one set for 
the first electron, and another set for the second electron. 
Using steepest descent methods, each set on its own can 
be used to calculate separate momentum yields of the 
form 
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for the first electron, and 



M('^(P2) 
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for the second electron (and we indeed do so in the fol- 
lowing section). However, the total RESI transition am- 
plitude Eq. ^ generally does not factorize in this way. 



M(Pi,P2)y^Af(i)(pi)M(2)(p2). 



(19) 



The reason is the time constraint in the original expres- 
sion ([T|) . In this original integral, as well as in fully classi- 
cal theories with real- valued trajectories, the times t", t' , 
and t are all real and can be ordered. However, the semi- 
classical pathways all have complex times. Therefore, the 
underlying issue of causality requires closer attention. In 
the remainder of this paper, we address this issue by ex- 
amining the key technical step in the derivation of semi- 
classical expressions, namely, the construction of steepest 
descent contours in the complex-time plane. In order to 
prepare this discussion, we set out in the following Sec. 
mil by providing semiclassical expressions for the partial 
amplitudes ^ and P^ . 



III. PARTIAL MOMENTUM-SPACE MAPS 

In this section we provide a detailed analysis of the 
partial one-electron transition probabilities |A/(^)(pi)|'^ 
and |Af(^)(p2)|^, as functions of the momentum compo- 
nents (j>i^2\\iPi,2Jl)- We evaluate the partial probabilities 
asymptotically via the steepest-descent method, which 
delivers expressions in terms of the saddles of the rapidly 
fluctuating integrals in Eqs. (IT71) and (|18p . respectively. 
The momentum-space maps obtained in this way will be 
useful for the construction of the correlated RESI tran- 
sition amplitude ^, Sees. IIVI and IV l in which causality 
must be taken into account. 

For M^^)(p2), the saddles are always well-separated 
in the momentum ranges of interest. Hence, it suffices 
to employ the standard saddle-point approximation, in 
which the saddles are treated individually and one per- 
forms a gaussian expansion around each of them. For 
the partial transition amplitude M'^'^^(pi) related to the 
first electron, however, the separation condition of sad- 
dles is not always fulfilled. In fact, at the borders of the 
classically allowed region in momentum space the saddles 
will approach each other closely. Hence, each pair of sad- 
dles must be treated collectively employing the uniform 
approximation discussed in [Sj. Beyond the boundary 
of the classically allowed region, one of the two saddles 
in a pair will lead to an exponentially increasing transi- 
tion amplitude and thus must be left out. This switching 
from two to one saddle near the boundary is known as 
a Stokes transition j25|, and occurs via a bifurcation of 
the steepest descent contour, which visits both saddles 
before the transition and only one saddle after the transi- 
tion. This happens when, for two saddles a and h defining 



a pair, Re[Si{pi,k,t'^X)] = Re[^i(pi,k, *;,,<;')]■ If, on 
the other hand, no classically allowed region exists, one 
must take the saddle that leads to exponentially decaying 
contributions throughout and the standard saddle-point 
approximation. (Similar Stokes transition will be crucial 
for the analysis of causality in Sec. II VI ) 



A. Monochromatic driving field 

Figure [1] shows the resulting partial transition proba- 
bilities |M(i)(pi)P and \M'-^'^ {p2)\'^ as functions of the 
perpendicular and parallel momentum components of 
each electron, for a monochromatic driving field. 

We start with the simpler case of the second elec- 
tron. For this electron, there are two saddles, whose 
start times, for vanishing parallel momenta, lie at ad- 
jacent field maxima separated by half a cycle of the driv- 
ing field. Each saddle on its own gives identical contri- 
butions to the probability |M*^^)(p2)|^. These individ- 
ual contributions are depicted in Fig. [IJc). Figure [Ijd) 
shows the rich interference pattern which arises from the 
quantum mechanical superposition of the two adjacent 
saddles. The individual contributions are maximal at 
{P2\\,P2±) = (0,0), for which the effective potential bar- 
rier is narrowest and tunneling most probable. These 
momentum components mark the center of the momen- 
tum map. The probability density quickly drops off with 
increasing momenta P2, i.e., as one moves away from 
the center of the momentum map. This is expected as 
the effective potential barrier through which the electron 
tunnels widens in this case (see discussion of Eq. p^ ). 

For the first electron, there exist two pairs of sad- 
dles. Each pair stems from a half cycle of the field 
and, for the parameter range employed in the figure, 
can be associated to a classically allowed region cen- 
tered at {p2\\,P2±) = (t2\/C^, 0) (see discussion of 
Eq. (iMl) ). Figures [TJa) and (b) show the resulting proba- 
bility |M^^)(pi)p, calculated either only using the saddle 
leading to exponentially decaying contributions outside 
the classical boundary [panel (a)] or the complete pair 
[panel (b)]. The boundary of the classical region is vis- 
ible as an outer, bright ring, which indicates the locus 
of the Stokes transition. If in panel (b) the second sad- 
dle would be dropped abruptly beyond the transition one 
would obtain a cusp. We eliminated this artifact by treat- 
ing the pair collectively using the uniform approximation 
[8| . The remaining fringes are caused by the interference 
between the two orbits of the pair. In practice, there 
is only interference between saddles separated by half a 
cycle around {pi\\,pi±) — (0,0), as they mostly populate 
different momentum regions. 



B. Few-cycle pulse 

We will now turn to the few-cycle pulse given by 
Eq. ([7]), with N = 5 and carrier envelope phase 4> = 0. 




FIG. 1: Partial REST transition probabilities for a monochro- 
matic driving field, as function of the momentum components 
Pi,2|| parallel and pi,2x perpendicular to the laser field polar- 
ization. Panels (a) and (b) depict the transition probability 
jAf'^'(pi)|^ for the first electron, Eq. ((TTl), while panels (c) 
and (d) give the transition probability |M'^'(pi)|'^ for the sec- 
ond electron, Eq. (|18p . Panel (a) shows the contribution of the 
dominant trajectories of the first electron (one per half cycle), 
while panel (b) shows its interference with the subdominant 
partner trajectories in the same half cycle, treated collectively 
in uniform approximation. The trajectory pairs from the two 
half cycles result in identical distributions for opposite paral- 
lel momenta; interference between these pairs of trajectories 
is negligible. Panel (c) shows the contribution of an individ- 
ual orbit of the second electron, while panel (d) shows its 
interference with second orbit in the same field cycle. The 
distributions have been normalized with regard to the largest 



value in each panel. Parameters are for Helium {El 



(9) 



0.97 



a.u., E^f' — 2 a.u. and E2 ~ 0.5 a.u.), and the driving 
field is monochromatic with intensity 7 = 3 x 10^*W/cm^ 
and frequency lj — 0.057 a.u. 



This pulse is shown in Fig. [5J together with the approxi- 
mate ionization and rescattering times for the first elec- 
tron, and the ionization times for the second electron 
[panels (a) and (b), respectively]. 

The first electron will leave close to an extremum of 
the field and return near the subsequent field crossing 
[Fig-E^a)]. The times indicated in the figure are associ- 
ated to the real parts of the solutions of the saddle-point 
equations pT|) - p5|) . for pi± = and pi[| = —A{t'). Each 
of the five arrows is actually associated to a pair of com- 
plex saddle points in the {t",t') plane. These pairs will 
be referred to as Pairs l(ei) to 5(ei). Upon recoUision, 
the returning electron will excite a second electron, which 
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FIG. 2: Schematic representation of the electric field E(f) and 
the corresponding vector potential A(t), for a few-cycle pulse 
of five cycles. The arrows in panel (a) indicate the approxi- 
mate times around which the first electron leaves, in case it 
returns at a crossing. The complex return and start times for 
the indicated pairs of orbits will have real parts in the vicinity 
of such times. The rectangles in part (b) mark the approx- 
imate ionization times for the second electron associated to 
the orbits utilized in this work (Orbits f to 5). The fields 
have been normalized to E{t)/Eo and A{t)/Ao, where Eo, Ao 
denote the field amplitudes. 



will then leave near the subsequent field extremum. The 
orbits related to these maxima are labeled Orbit 1(62) to 
5(62). The remaining orbits, which are not represented 
here, will play a negligible role, due to the fact that the 
ionization times for the first or second electron lie too 
close to the trailing edges of the pulse. This implies that 
the corresponding ionization probabilities will be very 
small. Below, we discuss momentum-space maps for the 
above-mentioned few-cycle driving pulse. 

We will start by analyzing the transition probabil- 
ity |M(^)(pi)|'^ for the first electron as a function of 
(PilhPii)? for each of the relevant pairs of orbits indi- 
cated in Fig. [21 As in the monochromatic case, a reason- 
able insight into the relevant momentum range is pro- 
vided by the orbit in each pair whose contribution de- 
cays exponentially outside the classically allowed region. 
In Fig. [31 we display the contributions of these orbits. As 
an overall feature, the center of such maps is no longer 




FIG. 3: Contributions from specific sets of orbits 

to the momentum-resolved RESI transition probabilities 
|M'^'(Pi)|^ of the first electron, Eq. ((TTJ), for the few-cycle 
pulse of Fig. [21 Panels (a), (b), (c) and (d) correspond to 
Pairs 2, 3, 4 and 5(ei), respectively (to eliminate distorting 
interference effects, only the dominant orbit, which remains 
physical beyond the classical boundary, has been taken into 
account). The pulse has peak intensity / = 2.5 x 10^*W/cm^, 
frequency cj = 0.057 a.u. and carrier-envelope phase <j) = Q; 
the other parameters are the same as in Fig. [T] 



located at (pi||,Pi±) = {±2^JTjp, 0), but varies from cycle 
to cycle. This is due to the fact that this estimate, albeit 
valid for a monochromatic field, no longer describes a 
field crossing for a few-cycle pulse. In fact, the exact po- 
sition of a field crossing will depend on the pulse envelope 
and on the carrier-envelope phase. Still, the momentum 
maps remain centered around vanishing transverse mo- 
menta piA_. This is expected, as the effective widening of 
the excitation energy for pi_L ^ that can be inferred 
from Eq. p^ occurs regardless of the shape of the exter- 
nal driving field. 

Apart from that, the magnitudes of the partial prob- 
abilities |M'^-'^)(pi)P will vary from cycle to cycle. 
Whether the contributions of a certain pair will be promi- 
nent, irrelevant or even vanishingly small will depend on 
several issues. If, for instance, the field amplitude is 
large for a specific cycle, the ionization probability for 
the first and second electrons are expected to be large 
as well. Thus, the contributions from orbits starting at 
such times are expected to prevail. Apart from that, if, 
for a particular set of orbits, the kinetic energy of the 
electron upon return is much larger than the energy re- 



quired to excite the second electron, according to Eq. p3)) 
there will be an extensive classically allowed momentum 
region. This will also play a role in making the contri- 
butions of a particular set of orbits prominent. In the 
specific case presented here, the contributions from Pair 
3(ei) and Pair 4(ei), depicted in Fig. [3l[b) and (c), re- 
spectively, are comparable, and at least three orders of 
magnitude larger than those from the other pairs. This 
is expected, as the ionization times related to both pairs 
are very close to the center of the pulse (see Fig. ^. 
Hence, there is a high tunneling probability for the elec- 
tron at these times. Further inspection, however, shows 
that the classically allowed region related to Pair 3(ei) is 
larger. This is related to the kinetic energy the electron 
exhibits upon return, which is highest for this pair. The 
contributions of the remaining pairs are less relevant, as 
the ionization times are closer to the trailing edge of the 
pulse. For instance, the contributions from Pair 5(ei), 
displayed in Fig. [3jd), are several orders of magnitude 
smaller than the other pairs. This is due to the fact that 
for this specific pair of orbits, the first electron does not 
return with enough kinetic energy to excite the second 
electron and still reach the detector; hence, rescattering 
is forbidden throughout. We have verified that this also 
happens for Pair l(ei). For that reason, we do not in- 
clude its contributions in the present figure. 




P2x/[Upl 



FIG. 4: Contributions of individual orbits to tlie momentum- 
resolved RESI transition probabilities |M"'(p2)|^ for the sec- 
ond electron, Eq. (|18|l . for the same few cycle pulse and atomic 
parameters as in Fig. (3] Panels (a), (b), (c) and (d) depict the 
contributions from the orbits 1, 2, 3 and 4(e2), respectively 
(see Fig. [U. 

We will now analyze the partial ionization probabil- 
ity of the second electron by computing |Af^^)(p2)P for 
Orbits l-4(e2) individually. Fig. 0] shows the outcome 



of these computations as functions of the parallel and 
perpendicular momenta {p2\\,P2±)- Similarly to what 
has been observed for the first electron, the centers of 
the momentum-space maps shift away from the posi- 
tion {p2\\,P2±) — (0,0). This shift is due to the lack of 
monochromaticity of the driving field, which introduces 
an asymmetry around the field extrema and leads to a 
slight bias towards either the positive or the negative par- 
allel momentum region. This bias is more pronounced for 
orbits whose ionization times approach the trailing edge 
of the pulse. For the specific pulse considered in this 
work, this can be observed in the contributions from Or- 
bits 3(62) and 4(62) [Fig.|4ljc) and (d), respectively]. Still, 
for all cases, the momentum-space maps remain centered 
at p2± = 0. This is a consequence of the fact that, effec- 
tively, the potential barrier through which the electron 
must tunnel is narrowest in this case, regardless of the 
shape of the field (see discussion of Eq. ([TSl) '). 

Furthermore, the closer the ionization time is to the 
pulse center, the larger the contributions from the corre- 
sponding orbit will be. For the specific pulse considered 
here, for instance, the contributions from Orbit 1(62) and 
Orbit 2(62), displayed in Figs. Ufa) and (b), respectively, 
are at least one order of magnitude larger than those of 
the remaining orbits. In the figure, we do not include the 
contributions from Orbit 5(e2), as these are at least four 
orders of magnitude smaller than those of the remaining 
orbits. 

As in the monochromatic-field case, the contributions 
of the above-stated sets of orbits must be added coher- 
ently. This will lead to interference maxima and min- 
ima, both in the momentum maps and in the electron- 
momentum distributions. In Fig. [5l we display the mo- 
mentum maps for the first and the second electron (left 
and right panels, respectively). The upper and lower 
panels, respectively, exhibit the contributions of the dom- 
inant sets of orbits, i.e.. Pairs 3(ei) and 4(ei), and Orbits 
1(62) and 2(62), and the overall partial momentum maps, 
in which all orbits have been included. 

Figs.ISKa) and (b) exhibit two circular regions for which 
the partial probability density is non- vanishing. These 
regions are roughly centered at (pi||,Pij.) = {±2^yUp, 0). 
Slight displacements away from these points are again 
related to the lack of monochromaticity of the field. The 
negative parallel momentum region is dominated by the 
contributions of Pair 4(ei). In fact, due to the large 
tunneling probability associated with it, it leads to the 
brightest spot in these panels. This pair interferes with 
Pair 2(ei), which is temporally displaced by a full cy- 
cle of the driving field; this interference is visible as the 
substructure in Fig. [SJb). The positive parallel momen- 
tum region is dominated by Pair 3(ei). In both panels, 
one may identify well-defined annular fringes, which are 
caused by the interference between the long and the short 
orbit of Pair 3(ei). These orbits are temporally close, and 
their contributions are several orders of magnitude larger 
than those of the remaining pairs l(ei) and 5(ei). 

The interference scenario is different for the second 
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electron. In this case, orbits located near different 
half cycles of the field lead to contributions in over- 
lapping momentum-space regions, roughly centered at 
{P2\\tP2±) = (0,0). For instance, in Fig. [Sljc), inclusion 
of the dominant Orbits 1(62) and 2(62) already leads to 
a rich interference pattern. Additional substructures ap- 
pear if the remaining orbits considered in Fig. 2] are taken 
into consideration, as shown in Fig. [5jd). 




FIG. 5: Partial transition probabilities for a few-cycle pulse, 
as function of the momentum components p,j|| and Pn±- 
Panels (a) and (b) depict the partial transition probability 
|M^^'(pi)P for the first electron, while panels (c) and (d) de- 
pict \M^ (P2)| for the second electron. These are obtained 
by coherently adding the contributions of orbits addressed in 
Figs. [3]and[3] (the field and atomic parameters are the same as 
in these figures). In panels (a) and (c), this sum is restricted 
to the dominant orbits of pairs 3(ei) and 4(ei), as well as 
1(62) and 2(62), respectively. The distributions have been 
normalized with regard to the maximum probability density 
in each panel. 



IV. CAUSALITY FOR COMPLEX TIMES 

In this section, we describe how the fact that ionization 
of the second electron must be subsequent to the rescat- 
tering of the first electron reflects itself in the complex 
time plane, and exemplify this for a monochromatic driv- 
ing field and a few-cycle pulse. This construction forms 
the basis of the calculation of correlated RESI momen- 
tum distributions, which are presented in Sec. |Vl 

Causality in the classical sense means that the first 
electron rescatters with the core prior to the second elec- 
tron being freed. If both the rescattering time t' of the 
first electron and the ionization time t of the second elec- 
tron were real, this condition would simply require t' < t. 



This condition is also embodied in the integration range 
of the SFA expression ([T|) for the transition amplitude. 
However, because tunneling is not a classical process, 
the solutions of the saddle-point equations pTjl - dTSl) are 
complex. Hence, in order to evaluate the integrals in 
terms of these orbits one must determine how causality 
manifests itself in the complex plane. 

This can be done by reinspecting the steps that lead 
to the saddle point approximation of the RESI transition 
amplitude ([T]). This amplitude can be rewritten as 



M(pi,p2) 



(20) 



= / dt" dt' d^kV^';'lv^'U'^''-P''''^'"''^M'-^'^ip2;t'), 



whereupon the time-conditioned amplitude 



M(2)(p2;i') 



diI^j^^)e*^^(P^'*) 



(21) 



of the second electron depends on the rescattering time 
t' via the lower integration limit and is part of an inte- 
grand, in contrast to the partial amplitude M'^^)(p2) = 
M(2)(p2;-cx)), Eq. dUl), evaluated in Sec.lHll 

Let us assume that we have integrated out k via a 
saddle point approximation (which is exact if one ne- 
glects the k dependence of the form factors) and denote 
the resulting effective action as 5(pi,p2,i",i',t). We 
then can deform the integration manifold over the three 
times t" , t' and t into the complex plane. The inte- 
grand can be analytically continued and does not pos- 
sess any poles. Therefore, the value of the integral does 
not change by this deformation. In the steepest-descent 
method ^20] , this deformation is carried out such that the 
quickly oscillating part e^'^'^P^'^^'*'* '* -' of the integrand 
changes into a smooth expression with maxima around 
the saddle points. To this end one associates to each sad- 
dle {t'^,t'^,ta) a sheet on which ReS{pi,p2,t,t' ,t") = 
ReS'(pi, p2, in, iJi, i"). (In more than one dimension, 
these sheets are not unique, but this does not affect any 
results ^j|-) Different sheets can meet in zeros of the in- 
tegrand, since there the phase is not well defined. Sheets 
have to be joined such that the total deformed manifold 
starts and ends at the original lower and upper integra- 
tion limits of the integral, respectively. Therefore, one 
also needs sheets connected to these boundaries, which 
are again constructed such that the integrand decays 
rapidly as one moves away from the integration limits. 
Following this prescription, not all saddle sheets will be 
part of the total deformed manifold, which means that 
only a restricted set of all saddles is contributing to the 
final expression. The saddle point approximation follows 
by expanding S around the maxima or integration limits, 
so that one obtains simple integrals for each saddle point 
on the total contour, as well as each boundary sheet. 
These individual contributions are accurate approxima- 
tions if the saddles and integration limits are not too 
close to each other; otherwise one needs to employ uni- 
form approximations. 



21]), we 



In order to apply this procedure to Eqs. ([20 
proceed in three interhnked steps. 



(I) Consider the time-conditioned amphtude 
M^'^\p2',t'). The saddles in the complex t plane 
are the same as for the partial amplitude (|18|) . and 
therefore determined by condition (llSp . However, the 
steepest descent contour now depends on the integration 
limit t' , as we have to construct a continuous contour in 
the complex t plane that links it to the upper integra- 
tion limit (at real -f-oo). The contour therefore starts 
with a segment of fixed Re[S'2(p2,i)] = R.e[S'2(p2,i')]- 
This boundary segment has then to be linked with 
contour segments passing through saddles ta, for which 
Re[S'2(p2,i)] = R.e[S'2(p2,ia)]- Compared to the calcu- 
lation of the unconditioned amplitude Af'^^''(p2), this 
has two effects: (i) There is an additional contribution 
from the boundary segment linked to i', which we will 
neglect as it is related to the electron-impact pathway 
[Ujt and (ii) only a subset of saddles ta contributing to 
-^ (P2) will lie on the contour conditioned by starting 
at t' . Besides this selection criterion, the individual 
saddle point contributions do not explicitly depend on 
t' . However, the number of relevant saddles changes 
at Stokes transitions, which for a given saddle ta occur 
when Re[S2{p2,t')] = Re[S'2(p2,ia)]- (This condition is 
not sufficient; whether there is indeed a Stokes transition 
can be verified by constructing the explicit integration 
contour, as carried out below.) 

(H) Next, we apply the saddle point approximation to 
the remaining integrals over t" and t'. We will focus at 
situations not too close to a Stokes transition. Within a 
range in t' , the contributions from the saddle points of 
M'^^)(p2;i') will then not depend on t' . Therefore, this 
term can be treated as approximately constant, and will 
not affect the saddle point conditions (ITT]) - (|13I) for t" and 

t'M- 

(HI) This in turn means that we now can substitute 
the lower integration limit in M(^)(p2; t'), which was con- 
sidered general so far, with the values of complex saddle 
points t'f,. 

As an upshot, we find that for each saddle point tra- 
jectory (t'^, t'f^) of the first electron, only a certain number 
of saddle points ta of the second electron will contribute. 
These are the ones that lie on the continuous steepest 
descend contour in the complex t plane that starts at 
the rescattering time i(, of the first electron. Since the 
saddle-point values ij, depend on the momentum pi of the 
first electron, but the remainder of the complex t contour 
(and in particular the location of the saddles ta) depends 
on the momentum p2 of the second electron, this gen- 
eralized causality requirement results in additional cor- 
relations for the momentum distributions, which will be 
quantified in the following Sec. [V] In the remainder of 
this section we focus on the explicit construction of the 
described steepest descent contours. 
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FIG. 6: Construction of steepest descent contours for the eval- 
uation of M^^'{p2;t'), Eq. H2ip . for monochromatic driving, 
with complex rescattering times t' corresponding to pi± — 0, 
Pi|| — 0. In panel (a), these rescattering times, obtained by 
solving the saddle point equations (|ll|) - (|13 p. are indicated by 
the star symbols. For orientation, the curves show the trace 
of times i' when piy is varied from — S-y/t/p to +51^1^ while 
Pi± = remains fixed. In panels (b)-(d), the circles indi- 
cate the ionization times t of the second electron for p2± — 0, 
while the parallel momentum takes the values P2|| =0 (panel 
b), P2\\ = 1.5^/Up (panel c), and p2\\ = S^/Up (panel d), 
respectively. These times are obtained by solving the saddle 
point equation (|15|) . The solid circles indicate saddles that are 
physically allowed (with positive Im5'2), while the open cir- 
cles are unphysical. The dark solid lines indicate the steepest 
descent contour through these saddles, while the dark dashed 
lines indicate steepest ascent contours, as well as steepest de- 
scent contour through unphysical saddles. The red (light) 
lines are steepest descent contours passing through the rescat- 
tering times i' (star symbols). These lines are all determined 
from the condition Re 52 = const. A full steepest descent 
contour must make use of solid lines in order to connect a 
lower integration limit t = t' (stars) to the upper integration 
limit (at t — -l-oo). The thick lines show the full contour for 
one selected starting value, indicated by the solid star. This 
contour only visits a selection of all physical saddles (dark 
dots). Between panels (c) and (d) a Stokes transition occurs, 
where an additional saddle becomes accessible. 



A. Steepest descent contours for monochromatic 
driving 

Steepest descent contours for monochromatic driving 
and different combinations of momenta are shown in Figs. 
[B][H1 In order to identify the starting point of such con- 
tours, we show in panel (a) of each figure the rescattering 
times t' of the first electron (star symbols) , corresponding 
to (Pii|,Pi±) = (0,0) in Fig.H ipn,pii_) = {2JU;,0) 



in Fig. [71 and {pi\\,pi±) 



, 0) in Fig. E (ac- 
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FIG. 7: Construction of steepest descent contours as in Fig. 
[Bl but for rescattering times corresponding to pi± — 0, 
Pi|| — 2^JTLp. Due to this change, in every cycle one pair of 
rescattering times moves closer to the real axis, corresponding 
to physically allowed rescattering, while the other pair moves 
away from the real axis. The initial values of the steepest de- 
scent contours in panels (b)-(d) are therefore shifted. For the 
rescattering time indicated by the solid star, again a Stokes 
transition occurs, but now a saddle is lost, not gained. 
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FIG. 8: Same as Fig. [6l but for rescattering times correspond- 
ing to pix = 0, Pi|| — —2^JUp. In comparison to Fig. [71 the 
rescattering times are shifted by a half cycle. The rescattering 
time indicated by the solid star now corresponds to a negli- 
gible contribution (its partner saddle is unphysical) and does 
not exhibit a Stokes transition. 



cording to Fig. [1] these values corresponds to the point 
niidvifay between the maxima in the momentum map, as 
well as the locations of these maxima). Panels (b)-(d) of 
Figs. in][5] show these values as integration limits in the 
t plane, along with the complex ionization times of the 
second electrons (dots). In all panels p2i_ = 0, while 
the parallel momentum takes the values P2\\ =0 (panel 
b), P2\\ = l-5y/Up (panel c), and P2\\ = ^\/Up (panel 
d), respectively [according to Fig. [U panel (b) therefore 
corresponds to the maximum in the momentum map of 
the second electron]. The solid lines in these plots show 
steepest descent segments from which a full contour must 
be constructed. These only pass saddles in the upper half 
plane (solid dots) , as the saddles in the lower half plane 
(open dots) have ImS'2 < 0. They will lead to expo- 
nentially increasing contributions and therefore are un- 
physical. However, for a given starting value not all the 
physical saddles are visited. This is illustrated by sample 
contours (thick lines) that start at a selected rescatter- 
ing time, indicated by the solid star, and only pass the 
saddles indicated by the dark dots, while the gray solid 
dots are left out. 

In this way, for any given saddle-point rescattering 
time t' we can define a clear boundary in the complex 
t plane that determines whether a saddle-point ioniza- 
tion time is allowed by causality or not. This boundary 
passes over saddles whenever a Stokes transition occurs. 
The thick sample contours display a Stokes transition in 
Figs, m and H but not in [51 

Inspecting these figures more generally, we see that to a 
good approximation, the causality requirement assumes 
the form Rei' < Rei. This is the case because Stokes 
transitions take place when a start time t' crosses over 
the steepest- ascenf segment of a physically allowed ion- 
ization time t. In the figures, these lines are seen to lead 
almost vertically from a physical saddle in the upper half 
of the complex-i plane to an unphysical mirror saddle 
in the lower half plane. We verified that this approxi- 
mate causality criterion remains valid over the range of 
momenta where the momentum maps in Fig. [Hare large. 



B. Steepest descent contours for a few-cycle pulse 

Figures [MTTI show the corresponding steepest descent 
contours for the few cycle pulse. Instead of a periodic 
repetition cycle-by-cycle, the rescattering times (stars) 
and ionization times (dots) are now modulated accord- 
ing to their position in the pulse (see Sec. IIII Bl and Fig. 
^. In the tails of the pulse, where the intensity is small, 
the rescattering and ionization times move away from 
the real axis. Compared to the monochromatic case, the 
steepest descent contours are therefore distorted, but the 
main features are still strikingly similar. In particular, 
we observe the same type of Strokes transitions, and the 
simplified causality criterion Ret' < Rei is clearly still 
a good approximation. However, from these figures we 
can still anticipate that the causality constraint in itself 
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FIG. 9: Construction of steepest descent contours as in Fig. 
[HI but for a few cycle pulse. Labels refer to the orbit pairs 
(first electron) and orbits (second electron) indicated in Fig. 
[21 Panel (a) shows the effect of the finite pulse length on the 
first electron. The stars again denote the rescattering times 
for pix ~ 0, pi|| = 0. In the tails of the pulse the driving 
is weak, resulting in rescattering times that are shifted away 
from the real axis. Therefore, some of these trajectories of 
the first electron become negligible. Panels (b)-(d) show that 
the same effect also occurs for the ionization times t of the 
second electron. The solid star indicates a rescattering time 
in the center of the pulse, which gives the largest semiclassi- 
cal contribution. The value of this time is close to that in the 
monochromatic case, and the associated steepest descent con- 
tours (thick lines) display the same type of Stokes transition 
as in that case (cf. Fig. [B]). 



plays a much more significant role when it comes to the 
calculation of actual transition probabilities. Large con- 
tributions to the momentum yield should stem from the 
saddle points close to the center of the pulse. Because of 
the causality constraint, however, it is not always possible 
to combine the most favorable rescattering and ionization 
times. In the figures, this is illustrated by the sample con- 
tours, which start at the rescattering time t' of the long 
orbit in pair 4(ei) (the most favorable time according to 
Fig. [3]). The most favorable ionization time for the sec- 
ond electron is that of orbit 2(62) (see Fig.H]), but these 
orbits generally cannot be combined because of causality. 
In the following section, we quantify the consequences in 
terms of the correlated momentum distributions of both 
electrons. 



FIG. 10; Construction of steepest descent contours for a 
few cycle pulse as in Fig. [9l but corresponding to pi± = 0, 
Pi|l = 2^/Up, for which the rescattering time indicated by the 
solid star moves closer to the real axis. In comparison to the 
monochromatic case in Fig. \7\ the associated steepest descent 
contours again display the same type of Stokes transition. 
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FIG. 11: Same as Fig. [9l but for rescattering times corre- 
sponding to pix ~ 0, Pi|| = —2y^7jp. As in the monochro- 
matic case (Fig. [8|, the rescattering time indicated by the 
solid star now corresponds to a negligible contribution and 
does not exhibit a Stokes transition. 
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V. MOMENTUM DISTRIBUTIONS 

In this section we compute two-electron momentum 
distributions as functions of the momentum components 
(Pi||jP2||) parallel to the laser field polarization, contrast- 
ing again the case of a monochromatic field to that of 
a few cycle pulse. For both cases, we calculate the mo- 
mentum distribution for resolved parallel and restricted 
ranges of the transverse momenta. If transverse momenta 
are fully integrated out, as it is done in many NSDI 
studies, quantum-interference effects get washed out and 
causality-related effects can no longer be identified. This 
also holds for cusps or further artifacts that may be 
present in saddle-point approximations, and which in- 
dicate their breakdown. Hence, transverse momentum 
integration would mask the very effects we intend to an- 
alyze. 

We consider the distributions 



^(Pl|h-P2||) 






Pi. 



P2_ 



,ax) 

'2„ j2 






d pij_d P2±|M(Pl,P2)+Pl ^ P2 

(22) 

which have been symmetrized with regard to electron 
exchange. This is necessary as the two electrons are in- 



(min) 



and 



distinguishable. In the above-stated equation, p^j^ 

Pn± (n = 1,2) denote the minimal and the maximal 
transverse momenta to be taken into account. 

For a monochromatic driving field, the causality- 
induced shielding of saddles does not have an observ- 
able effect on the momentum distributions. This is so 
because the saddles are repeating periodically cycle by 
cycle, and all have to be added coherently,. This even- 
tually results in a closely spaced train of delta func- 
tions which embody the Bohr resonance condition for 
multiple-photon absorption [26|. Smoothing over these 
delta functions, we find M(pi,p2) = M'-^'){pi)M'^^\p2) 
as if causality had been ignored. For a few-cycle pulse, 
however, the orbits in different cycles are not equiva- 
lent, and casuality has observable consequence. This 
is demonstrated in Fig. 1121 where we present the dis- 
tributions F{piH,p2i\) for a monochromatic field (upper 
panels) and a few-cycle pulse (middle panels ignoring 
causality, and lower panels respecting casuality), for in- 
creasing transverse momenta (left panels: 'low momenta' 
< Pi±,P2± < 0.2^,/U^; middle panels: 'medium mo- 



menta' 0.6 a 



<Pi 




<0.8, 



<Pi±,P2± < 1-2^ 



right panels pan- 



els: 'large momenta 

For the monochromatic driving field. Figs. fT2l a)-(c). 
all distributions adhere to the momentum constraints de- 
fined for RESI [2l|, i.e., they are cross-shaped and located 
along the axes p„|| — {n = 1, 2). The elongation of such 
distributions is determined by the saddle-point equation 
([HI) , which yields the momentum-space regions filled by 
the rescattering of the first electron, and the width of 
such distributions is given by the maximal and minimal 
momenta determined by Eq. (|16p . They are also sym- 
metric upon the refiectionp„|| -^ — p„||, which is a direct 



consequence of the time-reversal symmetry of the field. 
This means that, without symmetrization upon electron 
exchange, these distributions would be located along the 
horizontal axis pi|| = 0. The symmetrization pi -f-)- p2 
leads to the vertical axis of the crosses. The rather rich 
substructure present in the figure is associated to the 
quantum interference between the long and short orbit 
for the first electron, and between orbits displaced by 
half a cycle for the second electron. 

For a few-cycle pulse, in contrast, the electron mo- 
mentum distributions are very much affected by causal- 
ity. To isolate these effects, we start by analyzing the 
situation if causality could be neglected, displayed in 
Figs. [HKd), (e) and (f). These distributions are again 
obtained by assuming Af(pi,p2) = M'-^\pi)M'-^'){p2). 
The distributions are then asymmetric, which is a di- 
rect consequence of the asymmetries in the momentum 
maps in Fig. [5] The momentum maps are localized in 
the momentum region determined by the dominant or- 
bits, located close to the center of the pulse. For low 
transverse momenta, the contribution from Pair 4(ei) 
dominates in M'^^^(pi), which in Fig. [5ja,b) gives rise 



to the bright spot at {pi\\,Pii_) = (—2-^ 



Hence, 



one expects the distribution to be located in the region 



<P2\\ 



'Up along the negative pi|| half axis. 



Symmetrization upon electron exchange leads to the oc- 
cupation of the momentum-space region around the neg- 
ative P2\\ half axis, with —I^JTTp < pin < 2y^Up. 

For medium transverse momenta the contribution from 
Pair 3(ei) [corresponding in Fig.[5Ka,b) to the bright ring 
for pi|| > 0] becomes comparable to that of Pair 4(ei). 
Consequently, the two-electron distribution in Fig. H^ e) 
spreads out in the parallel momentum plane. Indeed, 
the contributions around the positive half axes p„|| > 
(n = 1, 2) are now comparable to those along the negative 
half axes. For large transverse momenta (Fig.[T2l(f)), Pair 
3(ei) dominates over 4(ei), and the distributions moves 
into the half positive axes. Apart from that, the results 
show that, the larger the transverse momentum range is, 
the more concentrated around p„j| = 2y/Up the electron 
momentum distributions are. 

Up to the present stage, however, causality has not 
been taken into account. For instance, the largest con- 
tributions to the partial momentum maps for the first 
and second electron [Pair 4(ei) and Orbit 2(62), respec- 
tively] , are not connected by causality, and therefore their 
combined contribution must be discarded. Indeed, the 
most favorable overall RESI pathway arises from com- 
bining Pair 3(ei) with Orbits 2(62) (if accessible) and 
3(62). The large contribution related to Pair 4(ei) is not 
sufficient to counter-act the lower tunneling probabilities 
related to Orbits 4(e2) and 6(62). Figs. [T2lg) to (i) il- 
lustrate the consequences. Since Pair 3(ei) results in a 
large yield at piu > 0, the electron momentum distri- 
butions are now mostly concentrated along the positive 
half axes Pn\\ (ji = 1,2). This situation persists over all 
transverse momentum ranges. 

Apart from the reshaping of the distributions, a note- 
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FIG. 12: Correlated RESI two-electron momentum distributions ^(pi||,P2||) as functions of the electron momentum components 
p„ll (n — 1,2) parallel to the laser-field polarization. The atomic parameters correspond to Helium (see Fig. [!}. Panels (a) 
to (c) have been computed for a monochromatic field of intensity / = 3 x 10^*W/cm and frequency a; = 0.057 a.u., while 
in panels (d) to (i) a five-cycle pulse of intensity / = 2.5 x 10^*W/cm , frequency lj = 0.057 a.u. and carrier-envelope phase 
(j> = has been used. The middle row of panels (d)-(f) show results for the few-cycle pulse if causality is discarded (these 
results are therefore unphysical), while the bottom panels (g)-(i) show the results for the few-cycle pulse including the causality 
restrictions. From the left to the right column of panels, we integrated over a transverse momentum range centered around 
low, medium, and large transverse momenta, respectively. 



"worthy feature in Fig. [12] is an overall decrease of a few 
orders of magnitude in comparison to the situation in 
"which causality has been disregarded. Taken altogether, 
these results sho"w that causality has a drastic effect on 
the two-electron momentum distributions in the RESI 
pathway. 



VI. CONCLUSIONS 

In this work we performed a detailed analysis of the 
recoUision with subsequent tunneling ionization (RESI) 
mechanism in laser-induced nonsequential ionization 
(NSDI), with emphasis on the implications of causality. 
Physically, RESI means that the first electron rescatters 
with the core at an instant t' and excites a second elec- 
tron, which tunnel ionizes at a subsequent time t. Causal- 
ity means that tunnel ionization of the second electron 
can only occur after the recollision of the first electron. 
Applying the saddle point approximation to the strong- 
field expressions, the rescattering and ionization times 
become complex (since tunneling does not have a clas- 



sical counterpart), and the notion of causality has to be 
generalized into the complex time plane. We have shown 
that the concepts of "before" and "after" translate into 
boundary conditions limiting the steepest-descent con- 
tours to be taken. These boundary conditions are given 
by the complex rescattering times of the first electron. 
Ionization times of the second electron are casually con- 
nected to a rescattering event if it lies on a steepest de- 
scent contour that connects the rescattering time to the 
distant future. In practice, this often translates into a 
simple rule where only the real parts of the complex times 
have to be compared. Deviations have been observed, 
however, for pairs of electron orbits associated with the 
trailing edge of the pulse, or for momentum regions in 
which the RESI yield is strongly suppressed. 

We illustrated the influence of causality and quantum- 
interference on momentum-resolved two-electron distri- 
butions. In order to isolated these effects, we compared 
distributions for a monochromatic driving field to those 
for a few-cycle pulse. While causality does not affect the 
electron momentum distributions obtained for monochro- 
matic fields, it significantly affects the results for a few- 
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cycle pulse. For these pulses, there exist several compet- 
ing sets of orbits, whose dominance influences the shapes 
of the electron momentum distributions, which generally 
are highly asymmetric and concentrated in specific mo- 
mentum regions. Causality puts constraints on the con- 
struction of two-electron trajectories, which drastically 
changes the momentum regions populated by the distri- 
butions. 

Causality may have consequences even for monochro- 
matic driving when it is combined with additional effects, 
such as bound-state depletion. Depending on the laser- 
field intensity and the binding energy of the excited state, 
such depletion can be considerable, and if this is the case 



the ionization time t lying closest in the future of the 
rescattering time t' is expected to contribute most. Such 
effects go beyond the present work, but may warrant fur- 
ther investigation. 
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The length-to velocity gauge transformation will intro- 
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However, the neglected contribution from the boundary 
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